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ABSTRACT 

'’A  comprehensive  theoretical  model  was  formulated  to  study  the  transient 
ignition  mechanisms  of  solid  propellants  under  rapidly  varying  pressure  con¬ 
ditions.  The  model  considers  a  two-dimensional  geometry  to  describe  the  ig¬ 
nition  of  a  heterogeneous  propellant.  The  rate  of  pressurization  has  been 
included  in  the  energy  equation  to  simulate  rocket  motor  ignition  conditions. 
The  model  does  not  make  any  a  priori  assumption  of  solid-phase,  heterogeneous, 
or  gas-phase  ignition  mechanism.  Chemical  reactions  are  considered  in  the 
gas  phase,  at  the  interface  as  well  as  in  the  subsurface.  The  most  up-to- 
date  chemical  kinetics  information  is  included  in  the  model.  A  generalized 
ignition  criterion  has  been  proposed  which  allows  ignition  to  occur  at  any 
site. ■ 

Development  of  the  theoretical  model  was  preceeded  by  an  extensive  review 
of  the  literature  on  solid-propellant  ignition  to  establish  the  state-of- 
the-art.  Various  ignition  theories,  experimental  measurements,  and  ignition 
criteria  were  critically  examined.  Major  technological  gaps  in  solid-propel¬ 
lant  ignition  study  were  identified,  and  areas  for  future  research  were 
recommended . 

A  study  was  conducted  to  investigate  flame-spreading  and  combustion  pro¬ 
cesses  in  narrow  propellant  cracks.  Experimental  results  show  that  when  the 
crack-gap  width  is  very  small  (~450  Up) ,  several  partial  closures  of  the  crack 
gap  occur  due  to  the  propellant  deformation  caused  by  strong  interaction  be¬ 
tween  combustion  and  solid  mechanics.  Results  from  the  theoretical  study 
of  the  effect  of  propellant  deformation  on  ignition  and  combustion  processes 
in  propellant  cracks  indicate  that  the  effect  of  propellant  deformation  is 
quite  important  for  small  gap  widths  and  high  chamber  pressure  conditions. 
Calculated  results  also  show  that  gap  closures  may  generate  substantial  local 
pressure  peaks. 


The  effects  of  pressurization  rate,  crack-gap  width,  crack  length,  and 
propellant  type  on  ignition  and  flame-spreading  processes  in  isolated  AP- 
based  solid  propellant  cracks  were  studied  experimentally.  It  was  observed 
that  hot  gases  precede  the  ignition  front,  and  that  the  maximum  pressure  in 
the  crack  cavity  may  significantly  exceed  the  pressure  in  the  chamber.  The 
results  of  parametric  study  indicate  that  as  the  pressurization  rate  or  burn¬ 
ing  rate  of  the  propellant  is  increased,  the  time  required  for  the  ignition 
front  to  reach  the  crack  tip  decreases,  and  the  maximum  velocity  of  the 
ignition  front  increases.  Maximum  pressure  in  the  crack  increases  with  an 
increase  in  burning  rate  or  crack  length,  but  decreases  with  an  increase  in 
gap  width.  _ _ _  ""I 
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I .  INTRODUCTION 


This  report  summarizes  the  progress  made  during  the  period  August  1,  1979 
to  July  31,  1980  under  the  project  entitled  "Transient  Ignition  Mechanisms  of 
Confined  Solid  Propellants  under  Rapid  Pressurization"  (Contract  No.  N00014- 
79-C-0762) . 

The  overall  objective  of  this  invetigation  is  to  achieve  a  better  under¬ 
standing  of  the  ignition  process  of  a  solid  propellant  under  rapidly  varying 
pressure  conditions.  These  fundamental  studies  are  expected  to  help  in  the  design 
of  solid  rocket  motors  and  in  reducing  the  hazards  caused  by  convective  burning 
and  deflagration-to-detonation  transition  processes.  The  report  covers  three 
major  areas:  1)  literature  survey;  2)  theoretical  work;  and  3)  experimental 
work.  These  topics  are  discussed  in  detail  in  Sections  II,  III  and  IV. 

The  following  publications  represent  a  part  of  the  work  performed  under 
the  contract. 

(1)  "Review  of  Solid  Propellant  Studies,"  AIAA  Paper  No.  80-1210,  AIAA/SAE/ASME 
16th  Joint  Propulsion  Conference,  June  30  -  July  2,  1980,  Hartford,  CT 

(by  A.  K.  Kulkami,  M.  Kumar,  and  K.  K.  Kuo). 

(2)  "Gas  Penetration,  Flame  Propagation,  and  Combustion  Processes  in  Solid  Pro¬ 
pellant  Cracks,"  AIAA  Paper  No.  80-1206,  AIAA/SAE/ASME  16th  Joint  Propulsion 
Conference,  June  30  -  July  2,  1980,  Hartford,  CT  (by  M.  Kumar,  S.  M.  Kovacic, 
and  K.  K.  Kuo) . 

(3)  "Combustion  -  Structural  Interaction  in  a  Viscoelastic  Material,"  Symposium 
on  Computational  Methods  in  Nonlinear  Structural  and  Solid  Mechanics,  October 
6-8,  1980,  Washington,  DC  (by  T.  Y.  Chang,  J.  P.  Chang,  M.  Kumar,  and  K.  K.  Kuo). 

(4)  "Effect  of  Propellant  Deformation  on  Igntion  and  Combustion  Processes  in 
Solid  Propellant  Cracks,"  to  be  presented  at  the  JANNAF  Propulsion  Systems 
Hazards  Subcommittee  Meeting,  October  27-31,  1980,  Naval  Postgraduate  School, 
Monterey,  CA  (by  M.  Kumar  and  K.  K.  Kuo). 
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II.  LITERATURE  SURVEY 

An  excellent  review  of  solid  propellant  ignition  research  was  conducted 

fourteen  years  ago  by  Price  et  al.^  (1966).  Numerous  ignition  studies  resulting 

in  significant  progress  have  been  made  in  the  interim.  To  fully  utilize  past 

2 

experience  and  knowledge,  an  extensive  literature  survey  was  made.  This  review 
was  presented  at  the  AIAA/SAE/ASME  16th  Propulsion  Conference  in  Hartford,  Con¬ 
necticut,  on  June  30  -  July  2,  1980;  the  abstract  of  this  paper  is  given  in  Appendix 
I.  Specific  objectives  of  the  review  were: 

(1)  to  establish  the  state-of-the-art  in  the  study  of  solid-propellant 
ignition; 

(2)  to  identify  technological  gaps  and  contradictions  in  the  existing 
literature; 

(3)  to  facilitate  development  of  a  comprehensive  ignition  theory; 
and 

(4)  to  incorporate  merits  of  past  experimental  techniques  into  the  current 
test  program. 

The  review  consists  of  an  easy-to-read  tabular  summary  of  recent  important 
works  in  the  area  of  solid  propellant  ignition,  giving  basic  assumptions,  theor¬ 
etical  formulation,  apparatus  description,  range  of  experimental  parameters, 
ignition  criteria,  important  results,  and  conclusions.  Critical  comments  on 
each  work  are  also  listed.  The  tabular  form  is  designed  to  facilitate  comparison 
between  various  studies  and  to  provide  a  quick  reference. 

In  general,  the  ignition  theories  can  be  classified  into  three  major  groups. 
Solid-phase  ignition  theories  emphasize  the  heat  deposition  rate  in  the  solid; 
the  ignition  criterion  is  usually  based  upon  the  attainment  of  a  critical  tem¬ 
perature  in  the  solid  phase.  Gas-phase  ignition  theories  consider  the  controlling 
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mechanism  to  be  the  chemical  reaction  between  vaporized  fuel  species  and  oxidizer 
gases;  the  ignition  criterion  is  usually  based  upon  the  rate  of  change  of  the 
gas-phase  temperature  distribution.  Heterogeneous  ignition  theories  propose 
the  controlling  mechanism  to  be  the  reaction  between  the  solid-phase  fuel  and 
ambient  oxidizer  at  the  interface;  the  ignition  criterion  is  usually  based  upon 
the  attainment  of  a  critical  surface  temperature  or  a  critical  rate  of  increase 
of  surface  temperature. 

Based  upon  study  of  over  100  investigations,  various  important  aspects  are 
discussed  in  the  review  in  order  to  present  an  overall  view  of  solid  propellant 
ignition  studies.  These  aspects  include  means  of  achieving  ignition,  description 
of  the  ignition  process,  governing  equations  and  interfacial  conditions  (together 
with  the  physical  meaning  of  each  term)  ,  effects  of  various  ambient  and  propellant 
parameters  on  ignition  delay,  difficulties  encountered  in  theoretical  and  experi¬ 
mental  investigations,  and  recommendations  for  future  research. 

As  a  result  of  this  review,  it  was  determined  that  some  of  the  major  tech¬ 
nological  gaps  in  the  theory  of  solid  propellant  ignition  are  a)  determination 
of  the  controlling  mechanism  at  various  stages  of  ignition,  b)  selection  of  a 
proper  ignition  criterion,  c)  detailed  specification  of  chemical  kinetics,  and 
d)  the  heterogeneous  nature  of  the  propellant.  Some  of  the  technological  gaps 
in  the  experimental  area  are  a)  realistic  simulation  of  an  actual  ignition  process 
in  the  laboratory,  b)  highly  transient  nature  of  the  ignition  process  and  the 
smallness  of  the  region  of  interest  (making  data  acquisition  very  difficult), 
c)  the  experimental  definition  of  the  instant  of  ignition,  and  d)  lack  of  broad 
data  base. 

It  is  concluded  that  although  the  foundation  for  solid  propellant  ignition 
research  is  laid,  it  is  as  yet  far  from  complete.  Additional  efforts  are  needed 
in  the  areas  of  a)  study  of  the  ignition  process  under  transient  conditions. 
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such  as  inclusion  of  pressurization  rate,  in  experiments  and  theory  (to  simu¬ 
late  actual  ignition  processes  in  rocket  motors),  b)  detailed  consideration  of 
chemical  kinetics  in  theoretical  models,  c)  consideration  of  heterogeneous  pro¬ 
pellants  using  2-D  or  3-D  models,  d)  use  of  advanced  instrumentation  in  conducting 
experiments  on  ignition,  and  e)  establishment  of  a  broader  data  base. 
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III.  THEORETICAL  WORK 

The  theoretical  work  was  divided  into  two  parts:  1)  development  of  a  compre¬ 
hensive  ignition  model;  and  2)  investigation  of  combustion/solid  mechanics  inter¬ 
action  in  narrow  cracks. 

3. 1  Development  of  a  Comprehensive  Ignition  Model 
3.1.1  Introduction 

Ignition  of  a  solid  propellant  is  a  complex  physicochemical  phenomenon  which 
involves  various  processes  such  as  heat  transfer,  fluid  mechanics,  phase  change, 
mass  diffusion  of  chemical  species,  and  chemical  kinetics.  Figure  1  shows  typical 

processes  involved  in  ignition  of  a  heterogeneous  propellant.  As  discussed  in 
2 

the  review  paper,  the  theoretical  modeling  of  the  ignition  process  is  complicated 
by  such  factors  as  the  heterogeneous  nature  of  the  propellant  (and  the  overall 
ignition  process),  specifications  of  chemical  kinetics,  determination  of  the 
controlling  mechanism  of  the  ignition  process,  and  selection  of  an  appropriate 
ignition  criterion.  In  existing  models,  a  number  of  simplifying  assumptions 
have  been  made  for  mathematical  tractability. 

Keeping  in  mind  the  difficulties  mentioned  above,  as  well  as  the  need  to 
explain  the  actual  ignition  event  adequately,  a  comprehensive  ignition  model, 
expected  to  be  the  most  complete  to  date,  was  formulated.  Highlights  of  this 
model  are: 

1.  2-D  geometry,  allowing  description  of  a  heterogeneous  propellant.  Thus 

3 

far,  only  Kumar  and  Hermance  have  used  a  2-D  ignition  model. 

2.  Inclusion  of  a  pressurization  (dp/dt)  term  in  the  governing  equation 
for  the  gases  surrouding  the  propellant  to  simulate  actual  rocket  motor 
ignition  conditions.  Up  to  now,  the  only  ignition  model  which  has 
considered  the  pressurization  effect  is  that  proposed  by  Kumar  and 

A 

Kuo  for  1-D  geometry. 


6 


SOLID 
PROPELLANT 


MELT 

LAYER 


GAS  PHASE 


External  heat  flux 
due  to  convection 
and  conduction 

Radiation  External 
absorption-^-radiant 
in  gas  heat 

phase  flux 


Exothermic?  Oxidizer 
gas-phase  j ^diffusion 
reaction  /  and  con¬ 


) 


vection 


Exothermic/ 


surface 

reaction 


Phase 
transition 


Variations  of 
ambient  conditions 
(e.g.,  pressuri¬ 
zation) 


Fig.  1  Physical  Processes  Involved  in  Solid 
Propellant  Ignition 


7 


3.  No  a  priori  assumption  of  a  solid-phase,  heterogeneous,  or  gas-phase 
reaction  mechanism  was  made  to  control  the  ignition  event.  Except 
for  the  unified  ignition  theory  of  Bradley,^  for  which  very  limited 
results  were  obtained,  this  approach  is  unlike  that  of  any  other  pre¬ 
viously  proposed  model.  The  current  model  allows  chemical  reactions 
in  all  regions,  including  the  gas  phase,  interfaces,  and  subsurface. 

4.  A  flexible  "multiple"  ignition  criterion  has  been  proposed  to  allow 
ignition  to  occur  at  any  site,  i.e.  ,  in  the  gas  phase,  at  the  interface, 
or  in  the  solid  phase. 

5.  Detailed  chemical  kinetics  information  was  considered  in  development 
of  the  model. 

3.1.2  Physical  Description  of  the  Ignition  Model 

To  model  the  ingition  of  a  heterogeneous  propellant,  one  must  first  examine 
the  nature  of  the  propellant.  Figure  2  shows  typical  side  cross-sectional  and 
top  views  of  a  composite  propellant.  Oxidizer  particles  of  random  size  (within 
a  certain  range)  are  distributed  in  a  fuel  binder.  For  a  mathematical  simulation 
of  this  type  of  structure,  a  statistical  size  distribution  of  oxidizer  particles 
is  estimated.  From  this  calculation,  the  radius  R^  of  the  oxidizer  particle 
and  the  radius  R£  of  the  fuel  binder  for  a  typical  element  of  the  propellant 
are  determined. 

A  modified  cross-sectional  view  of  a  statistically  averaged  element,  suitable 
for  mathematical  modeling,  is  shown  in  Figure  3.  A  cylindrical  coordinate  system 
(r,z)  is  selected  in  which  initially  the  r  axis  runs  along  the  surface,  starting 
from  the  center  of  the  oxidizer  particle  (which  is  assumed  to  be  hemispherical 
at  t**0) .  The  axis  is  perpendicular  to  the  initially  planar  interface  and  positive 
in  the  surrounding  gas  phase.  An  external  radiation  flux,  incident  upon  the 
interface  (see  also  Figure  1),  is  partly  absorbed  and  partly  reflected  at  the 
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PHYSICAL  MODEL 


A  Typical  Element 


Fig.  2  Schematic  Diagram  of  a  Heterogeneous  Solid  Propellant 
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Fig.  3  Modified  Cross-9ectional  View  of  a  Statistically 
Averaged  Element  Considered  in  the  Ignition  Model 
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interface,  and  partly  absorbed  in  depth  in  the  solid  phase.  Some  of  the  in-depth 
absorption  of  radiation  can  cause  photochemical  decomposition  and  some  radiation 
is  emitted  by  the  surface  to  the  surroundings.  Becuase  of  external  heat  deposition 
(by  radiation  or  by  hot  gases  in  the  surroudings) ,  the  fuel  and  oxidizer  begin 
to  gasify  when  the  gasification  temperature  of  each  individual  component  is  reached. 
The  gasified  fuel  and  oxidizer  diffuse  in  the  surroundings;  the  ambient  oxidizer 
may  diffuse  toward  the  interface  and  eventually,  following  the  mixing  of  these 
gases,  there  is  an  exothermic  reaction  in  the  gas  phase.  Simultaneously,  the  sub¬ 
surface  and  heterogeneous  reactions  may  continue.  Heat  evolved  and  absorbed 
from  all  of  these  reactions  is  redistributed  by  conduction,  convection,  radiation, 
and  mass  diffusion.  Eventually,  thermal  runaway  is  reached  at  any  of  the  solid, 
gas,  or  interface  sites.  The  ignition  criterion  then  determines  the  location  and 
instant  of  ignition.  The  above  process  is  described  by  the  following  mathematical 
formulation. 

3.1.3  Governing  Equations 

Equations  (1)  through  (4)  represent  the  gas-phase  continuity,  energy  and 
species  equations  for  oxidizer  and  fuel  gases. 


Continuity  equation: 
dP°  \  S(PgVz} 


3t  3  z 


Energy  equation: 


small 


cp  il+pcv  il-iP.I  II  +  J_  ,v  2L  , 

PPg  3t  Pg  p  z  3z  3t  r  JrAkgr  3r5  +  3z  (kg  3z)  +  qg 


Oxidizer  species  equation: 


(1) 


(2) 


3Y„ 


°8  '5T  +  V.  ~£r  -  r  A  +  “ox 


(3) 


i 
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Fuel  species  equation: 


3t 


3Yf  l  a 

pgVz  TT  =  7  3r  (rPpg  -3T> 


(4) 


The  inhomogeneous  terms  and  will  be  discussed  later. 

Since  the  region  of  interest  in  the  gas  phase  during  the  ignition  process 
is  very  small  (in  the  order  of  1  mm) ,  local  pressure  in  this  region  is  considered 
to  be  spatially  uniform.  However,  the  pressure  is  allowed  to  vary  with  respect 
to  time.  The  pressure-time  relationship  is  assumed  to  be  known  either  from  a 
measured  P-t  trace  or  a  prescribed  pressure-time  information  simulating  any  ig¬ 
nition  process.  Therefore,  the  gas-phase  momentum  equations  are  reduced  to 


P  -  P(t) 


(5) 


Becuase  pressure  is  usually  not  extremely  high  during  the  ignition  interval, 
the  perfect  gas  low  can  be  used  as  the  equation  of  state, 


P 


g 


_P 

RT 


(6) 


Becuase  of  the  heterogeneous  nature  of  the  propellants,  separate  energy 
equations  are  written  for  oxidizer  and  fuel  binder  in  the  solid  phase,  as  shown 
in  Equations  (7)  and  (8).  The  two  equations  govern  the  variation  of  temperature 
distribution  in  different  regions  with  different  thermal  properties. 


For  oxidizer: 


For  fuel  binder: 


where 


30x,sC0x,s  -57  ■  [k 


32T 


+  k„ 


3  T  .  k0x,s  3T ,  . 

+  — T1-  rrl  +  qr 


3t  Ox, s  ^2  Ox, s  ^r2  r  3rJ  40x,i 


0  „  il.rk  !?Z+k  i!l +  ^Ls  3T,  . 

PF , sc F , s  3t  [kF,s  3z2  +  kF,s  ^2  +  r  3r]  qF,s 


q"'  -  q"‘  +  q'»  +  a’" 

Ox,s  Ox, radiation  ^Ox, photochemical  Ox, pyrolysis 


(7) 

(8) 

(9) 
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-B0x  2 

q0x, radiation  T0xG^0xe  “  ^Qx^z 


(10) 


*Ox,  photochemical 


T0x  ^z^0x,pc^ri0x^d^ 


(following  Ohlemiller  and  Sumraerfield,  1968) 
q0x, pyrolysis  =  “  A0x,py  ^Ox,py  exp^-E0x,py^RuT) 


(ID 


(12) 


(following  Ohlemiller  and  Sununerf ield®  1968) 


3.1.4  Initial  and  Boundary  Conditions 

The  above  governing  equations  are  subjected  to  a  number  of  initial  and  boundary 
conditions.  The  following  table  summarizes  the  initial  and  boundary  conditions 
needed  for  the  solution  of  gas-phase  equations. 


Table  1  Initial  and  Boundary  Conditions  for  Gas  Phase 


If 

Equation 

■  - 

Major  Dependent 
Variable 

Independent  Variables 
for  BC  and  IC 

Subtotal  of 
IC  and  BC 

1 

Continuity 

V 

z 

t,z 

2 

2 

Energy 

T 

t,r,r,z,z 

5 

3 

Species:  Ox. 

Y0x 

t,r,r,z 

4 

4 

Species:  Fuel 

yf 

t,r,r,z 

4 

TOTAL  (IC  and  BC) 

15 

12 


13 


have 


14 


and 


=  Tl 

Z„  +  2_  - 

F-g  F-g 


(25) 


ZOx-g+ 


"Ox-g 


The  energy  flux  balance  across  the  oxidizer-gas  interface,  i.e.  z  =  z 


is  given  by 


Ox-g’ 


Ox,s  3z 


zOx-g 


-  k  ^ 
g  32 


+•  a 


(conductive 
heat  flux 
into  the  oxi¬ 
dizer) 


Ox-g 

(conductive 
heat  flux 
from  gas 
phase) 


Ox-gG  eOxEb 


'Ox 


(radiation  (radiation 
absorption  loss  from 
at  the  the  surface) 

surface) 


+  r  P  T(c.  -  c  )  + 
b0x  °X,S  °x,s  p’8 

(net  convective  heat  flux 
into  the  interface) 


MOx-g 

(net  heat  generation 
at  the  surface) 


(26) 


and  that  across  the  fuel-gas  interface,  i.e.  z  *•  z„  ,  is  given  by 

r-g 


v  il 

r,s  3z 


lZF-g+ 


k  & 
g  3z 


'F-g 


(conductive  heat  (conductive 
flux  into  the  heat  flux 


fuel) 


from  the  gas  at  the 


phase) 


+  aF-gG  -  GFEb 


(radiation  (radiation  loss 
absorption  from  the  surface) 


surface) 


+  r  p  T(c  -  c  )  +  q" 

b  F,s  F,s  p,g  MF-g 


(net  convective  heat 
flux  into  the  interface) 


(net  heat  generation 
at  the  surface) 


(27) 
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2  .  2 


At  the  oxidizer-fuel  interface, £=  /r  +  z  =  ,  the  energy  flux  balance  is 


-  K.  ^ 
^,s  35 


k  3T  . 

Ox,s  35  -  qOx-F 

R1 


(conductive  heat  (conductive  (net  heat  generation 

flux  into  the  heat  flux  at  the  surface) 

fuel)  from  the  oxi¬ 

dizer) 


+  a0x-FT0xG  e 


(radiation  absorption 
at  the  interface) 


2  ,  2 


At  the  5  =  /r  +  z  =  interface  ,  the  temperature  continuity  is 


The  mass  continuity  across  the  oxidizer-gas  interface  relates  the  local 
gas  velocity  to  the  regression  rate  of  the  oxidizer  by  the  following  relationship. 


0  v 

g  2  + 

ZOx-g 


POxrb„ 


similarly  for  fuel  we  have 


PgVz  +  "  PF\  • 

ZF-g  F 

3.1.5  Chemical  Kinetics  Information 

In  order  to  close  the  model  proposed  above,  the  following  terms  must  be 

expressed  in  terms  of  the  dependent  and  independent  variables  r.  ,  r,  , 

Ox  F 

qg"’  ^Ox  ’  “p"»  qOx-F’  qF-g’  “ox*  and  “f'-  F°r  thi®  PurP°se*  ^nowledge 
of  detailed  chemical  reactions  and  kinetic  data  is  necessary.  A  survey  of  kinetic 
information  was  therefore  conducted  to  obtain  thermal  decompositon  data  (rate 


i 

t 
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constants,  physical  constants  for  melting  and  ablation,  etc.),  and  to  determine 
the  key  steps  in  a  chemical  reaction  mechanism,  a  process  which  is  usually 
extremely  complicated.  In  particular,  kinetic  information  was  sought  for  AP, 
PBAA,  HMX,  RDX,  and  NC/NG;  a  summary  of  the  information  obtained  for  each  of 
the  ingredients  follows. 
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(a)  AP  (Ammonium  Perchlorate) 

AP  is  the  most  widely  used  oxidizer  in  composite  solid  propellants.  Most 

of  the  following  information  was  obtained  from  Guirao  and  Williams;^  however, 

8—11 

many  of  the  values  were  collected  from  and  confirmed  by  other  sources.  Briefly, 

the  overall  burning  of  AP  involves  first  the  dissociative  sublimation  into  ammonia 
and  perchloric  acid,  which  begins  around  473K  and  is  an  endothermic  process. 

The  perchloric  acid  then  decomposes  in  the  gas  phase  at  higher  temperatures  (above 
523K) ,  and  reacts  with  ammonia  exothermically.  The  products  of  this  reaction 
are  oxidizer-rich  gases  which  react  with  fuel  gases  vaporized  from  fuel  binders 
of  a  composite  propellant.  For  temperatures  between  473  and  573K,  the  thermal 

g 

decomposition  of  AP  can  be  represented  by 

4NH4C104  -*■  2  Cl2  +  3  02  +  8  H20  +  2  N20  (32) 

The  reaction  for  decomposition  above  623K  is 

80  NH.C10.  -*■  20  Cl.  +  16  N„0  +  20  N0C1  +  8  HC10.  +  12  HC1 
4  4  2  2  4 

+  75  H20  +  14  N2  +  51  02  (33) 

These  reactions  are  exothermic,  and  are  believed  to  be  preceded  by  dissociative 
7 ,9 

sublimation  ’  of  AP  as  shown  below. 


E  -30  kcal/mole 

NH.  clO.  -  ,  - - - - - - ► 

4  4,  .  (dissociative 

3  sublimation) 


NH-  +  HC10 . 
3Cg>  4<g) 


+  58 

mole 


(34) 


where 


+  HC10. 
4 


Oxidizer-Rich 

Product 

Gases 


-  772  kcal/kg  (35) 


d[HC104) 


dt -  k2  Ihcio4][nh3] 

E 


k-  -  A  T  exp  (- 

®2 


82 

R  T 
u 


) 


“  15.47  kcal/mole,  A  ■  2.03  x  10**  cm^/mole-s-k 

«2 


(36) 

(37) 


2 
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(b)  PBAA  (Polybutadiene  Acrylic  Acid) 

PBAA  is  a  fuel  binder  used  with  AP  in  many  composite  propellants.  PBAA 
is  usually  cured  with  EPON.  The  fuel  binder  made  with  86.6%  PBAA  and  14.4%  EPON828 

9 

is  represented  by  the  formula  CHj  $25° 0  03873'  thenna^  decomposition  of 

1 2 

PBAA  begins  around  610K  and  continues  until  770K.  According  to  Varney  and 
12 

Strahle,  the  decomposition  can  be  written  as 


PBAA 


Fuel-Rich  +  1950  kcal/kg 

[Product  Gases 


(38) 


d  In 


^BAA 


dt 


-  A  exp 


r. !_] 

I  R  T 

U  11  -i 


(39) 


where 


E 

A 


34  kcal/mole 
2  x  107  s"1 
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(c)  HMX  (Cyclotetramethylenetetranitramine) 

The  principal  products  of  HMX  decomposition  are  formaldehyde  (HCHO) ,  nitrogen 
oxide  (N20),  hydrogen  cyanide  (HCN) ,  and  nitrogen  dioxide  (NO^.13  The  currently 
understood  chemical  kinetic  processes  between  these  principal  products  are  given 
below. 


HCHO 

+  no2 

(fast) 

Products 

(AO) 

HCHO 

+  N2° 

(slow) 

Products 

(41) 

HCN 

+  no2 

(slow) 

Products 

(42) 

HCN 

+  N20 

7 

Products 

(43) 

Corresponding  rates  of  depletion  of  NO,,  and  N20  of  the  above  from  reactions  are 
given  below. 


d(N02) 

dt~ 

d(N20) 

dt 

d(N02) 

dt 

d(N20) 


1.26  x  1013  exp(26700/R^T)  (Ar) "4<3 (HCHO) ’ 56 (N02) ’ 90  mole/cm3  sec 
8.41  x  1011  exp(-27400/RuT)CAr)"'14(HCHO),64(N20)1,17  mole/cm3  sec 
4.9  x  10U  exp(-43300/RuT)(Ar)1'03(HCN)~'03(NO2) ,65  mole/cm3  sec 


dt 


*  not  known 


(d)  RDX  (Cyclotrimethylenetrinitramine) 

14 

According  to  Ben  Reuven  et  al.,  the  decomposition  of  RDX  can  be  represented 
by  the  following  reactions. 


RDX  |  N2  +  N20  +  N02  +  3  HCHO 

(gas  or  liquid) 

HCHO  +  N02  ^l0Wi  NO  +  y  C02  +  2  CO  +  y  H20 

HCHO  +  N20  — £*■  Products  (too  slow,  therefore,  neglected) 

(liquid)  «  3  x  1018  exp(-47500/RuT),  sec"1 
k^gas)  =  lO15,5  exp(-41500/RuT),  sec"1 
k2  =  101^exp(-19000/RuT) ,  cm3/gmol-sec,  T  >  430  K 

k2  =  1013'1  exp(-26700/RuT),  cm3/gmol-sec,  970K  <  T  <  1470K 


(44) 


(45) 


(46) 


(47) 


(48) 
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(e)  Double-Base  Propellants 

According  to  Kubota  et  al.,^  the  chemical  kinetic  processes  in  decomposition 
and  the  reaction  of  double-base  propellants  can  be  represented  by  the  following 
mechanism. 


(Condensed 

Phase) 


(Luminous 

(  urface)  •  (Fizz  Zone)  (Dark  Zone)  Flame  Zone) 


Propellant  ^-ac-tion€ 


Xl+Yl 


X1  +  V 


Reaction  Reaction 

(D  © 

(Gasification) 


x2  +  Y2 


X2  +V2 


Solid  Carbon 


Path  1  only 
Paths  1  and  2 

£l’e2 

A 

X,Y 

Z 


Eg  *  5  x  10  cal/mole 


Noncatalyzed  reaction  scheme 
Catalyzed  reaction  scheme 

Mass  fractions  of  the  propellant,  e^  +  =  1 

Products  after  decomposition 

Chemical  species  in  gas  phase 

Pre-exponential  factor 

6  3 

Za  a  5  x  10  cm  /g-sec 


Ev  ■  17  x  10  cal/mole 

1  3 

E,  *  17  x  10  cal/mole 

2  3 

E  *  17  x  10  cal/mole 

1  3 
E  *  7.7  x  10  cal/mole 


Z,  5  5  x  10  cm/sec  (?) 

1  5 

Z.  ”  5  x  10  cm/sec 

2 

Z  *  1  x  10  cm^/g-sec 

1  9  3 

Z  *  1  x  10  cm  ,  g-sec 

2 
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In  the  above  summary  of  kinetic  information,  some  of  the  detailed  inter¬ 
mediate  and  chain  reactions  are  not  given.  One  reason  for  omitting  them  is 
that  these  detailed  chain  reactions  are  usually  proposed  reaction  mechanisms 
which  have  not  been  verified  and  therefore  may  be  unreliable.  It  was  apparent 
from  this  survey  that  kinetic  information  in  general  is  incomplete  and  scattered. 
Because  of  the  inadequacy  of  relevant  kinetic  data  and  the  mathematical  intricacy 
involved  in  consideration  of  detailed  kinetics,  it  is  prudent  at  this  time  to 
consider  a  single-step,  global  reaction  between  the  oxidizer-rich  and  fuel-rich 
gases  for  determining  gas-phase  heat  release.  However,  realistic  rate  equations 
for  the  decomposition  of  oxidizer  particle  and  fuel  binder  are  taken  into  con¬ 
sideration  in  the  model. 

AP/PBAA  composite  propellant  is  considered,  in  the  solution  of  the  proposed 
ignition  model  for  the  following  reasons:  1)  a  number  of  experiments  have  been 
performed  and  are  in  progress  on  this  propellant;  and  2)  reliable  kinetic  infor¬ 
mation  is  available  for  this  propellant.  Numerical  solution  which  uses  AP/PBAA 
propellant  data  is  in  progress. 

3.2  Investigation  of  Combustion/Solid  Mechanics  Interaction  in  Narrow  Cracks 

Of  the  several  interdependent  gas  dynamic  and  solid  mechanics  processes 
that  may  influence  combustion  behavior  in  a  propellant  crack,  the  coupling  between 
the  pressure  along  the  crack  and  the  deformation  was  considered  to  be  the  most 
important.  The  analysis  considers  deformations  caused  by  both  pressure  loading 
and  burning.  Other  coupled  processes  such  as  thermal  stresses  due  to  the  temperature 
gradient  in  the  solid,  variation  in  burning  rate  due  to  compression  of  the  pro¬ 
pellant,  etc.,  were  considered  to  be  of  secondary  importance  and  were  ignored 
in  the  present  analysis. 
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The  general  layout  of  the  computer  program  is  shown  in  Figure  4.  The  theor¬ 
etical  coupling  considered  between  the  convective  burning  and  structural  analysis 
is  simplified,  as  discussed  previously,  to  the  transfer  of  pressure  loading  and 
crack  deformation.  Pressure  is  calculated  by  the  gas  dynamic  portion  of  the 
computer  program  at  each  nodal  point  location  on  a  one-dimensional  grid  along 
the  length  of  the  crack.  The  convective  burning  analysis  of  the  crack  combustion 
incorporates  the  crack  geometry  variation  caused  by  both  mechanical  deformation  and 
mass  loss  through  gasification  of  the  propellant  at  the  surface.  Once  the  gas- 
phases  equations  are  solved ,  and  pressures  and  burning  rate  along  the  crack  are 
calculated  for  a  particular  time  t,  program  control  is  transferred  to  the  nonlinear 
finite-element  analysis  program  (NFAP)  portion  of  the  combined  code.  NFAP  re¬ 
formulates  the  geometry  because  of  the  material  loss,  and  updates  the  stiffness 
matrix  for  the  new  time  step.  Surface  elements  of  the  finite-element  mesh  are 
loaded  with  the  pressure  obtained  through  the  gas-phase  equations.  The  propellant 
deformation  is  calculated  in  NFAP,  using  a  static  analysis  at  time  t.  A  plane- 
strain  analysis  is  used  to  solve  general  stress-strain  equations;  this  is  congruous 
to  the  experimental  test  configuration.  The  transient  nature  of  the  pressure 
loading  is  considered  by  using  a  static  analysis  at  incremental  time  steps. 

Figure  5  is  a  diagram  of  the  finite-element  grids  used  for  structure  analysis 
and  the  finite-difference  nodes  used  for  gas-phase  solutions.  The  configuration 
is  compatible  with  the  geometry  considered  in  the  theoretical  model,  as  well 
as  with  that  used  in  the  experiments.  Two-dimensional,  eight-node,  isoparametric, 
quadrilateral  elements  were  used  to  model  a  sample  crack.  Because  of  symmetry, 
only  half  of  the  crack  needs  to  be  considered.  Of  the  eight  nodes  In  an  element, 
four  are  located  at  the  corners,  and  one  at  the  midpoint  of  each  side.  Except 
near  the  tip  region,  the  finite-difference  nodes  used  for  the  gas-phase  solution 
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are  in  a  one-to-one  correspondence  with  the  finite-element  nodes  which  form  the 
propellant  surface  in  the  two-dimensional  structural  analysis  mesh.  This  is 
helpful  in  reducing  the  computational  effort  in  transferring  the  values  from 
finite-difference  nodes  to  finite- element  grids,  and  vice  versa. 

Figures  6  and  7  show  calculated  pressure  distributions  at  various  times 
for  cases  in  which  propellant  deformation  becomes  very  important.  For  results 
shown  in  Figure  6  (Case  1),  the  chamber  pressurization  rate  was  3.8  x  10^  atm/s 
and  the  gap  width  was  0.51  nan.  Only  those  curves  of  interest  are  shown,  and 
for  ease  of  explanation,  they  are  not  super-imposed.  Between  t  =  0.25  and  .275  ms, 
the  crack  gap  near  the  entrance  is  partially  closed,  causing  a  local  peak  near 
the  entrance;  however,  pressure  at  the  crack  entrance  is  identical  to  that  in 
the  chamber.  At  t  *  0.3  ms,  high  pressure  in  the  chamber  and  at  the  crack  entrance 
cause  the  closure  to  move  to  x/L  *  0.1,  resulting  in  a  local  pressure  peak.  The 
localized  high  pressure  at  x/L  *  0.1  causes  the  propellant  to  deform;  the  region 
of  gap  closure  moves  to  x/L  =0.2,  resulting  in  a  substantial  increase  of  local 
pressure.  The  subsequent  closure  results  in  a  zero  gap  width  and  causes  the 
gas-dynamic  solution  to  explode.  This  predicted  movement  of  the  partial  gap 
closure  is  very  similar  to  that  observed  experimentally  (see  Figure  11). 

Case  2,  shown  in  Figure  7,  was  computed  for  a  higher  pressurization  rate 
of  5  x  10^  atm/s  and  a  wider  initial  gap  width  of  .89  mm.  At  0.45  ms  (not  shown 
in  the  Figure),  gap  closure  occurs  near  the  crack  entrance,  resulting  in  localized 
high  pressure  near  the  entrance.  Again,  both  the  high  pressure  near  the  entrance 
in  the  cavity  and  that  acting  on  the  side  walls  exposed  to  the  chamber  cause 
the  gap  closure  to  move  downstream,  resulting  in  a  pressure  peak  near  x/L  =  0.17 
(at  .475  ms).  At  0.5  ms,  two  regions  of  gap  closures  are  obtained,  resulting 
in  pressure  peaks  at  x/L  =  0.1  and  0.27,  respectively.  As  time  passes,  the 


Fig.  6  Calculated  Pressure  Distribution  at  Various 
Times  during  Crack-Gap  Closure  (Case  1) 
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pressure  peaks  move  downstream;  some  disappear  as  a  result  of  propellant  defor¬ 
mation.  At  0.525  ms,  another  gap  closure  develops  near  the  entrance.  This  process 
continues  as  time  progresses.  Once  again,  these  results  explain  the  behavior 
of  the  gap  closures  observed  experimentally. 

Results  obtained  here  indicate  that  closure  of  the  crack  gap,  which  may 
initially  occur  at  the  crack  entrance  because  of  small  gap  widths  and  high  chamber 
pressures,  may  propagate  along  the  crack  and  result  in  local  pressure  peaks. 

The  gap  closure  is  not  observed  for  cracks  with  large  gap  widths  or  at  low  chamber 
pressures  because  propellant  deformation  under  these  conditions  is  small  and  does 
not  affect  the  ignition  process  substantially.  Detailed  governing  equations 
and  verification  for  the  structure  mechanics  part  is  given  in  Ref.  16.  A  paper, 
which  includes  a  more  detailed  analysis  of  the  combined  program,  has  been  accepted 
for  the  JANNAF  Propulsion  Systems  Hazards  Subcommittee  Meeting  (Ref.  17).  Abstracts 
of  these  papers  are  given  in  Appendices  II  and  III,  respectively. 
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IV.  EXPERIMENTAL  WORK 

4 . 1  Development  of  a  New  Igniter  System 

A  new  igniter  system  has  been  designed  and  fabricated  to  obtain  higher  pres¬ 
surization  rates  (>10^  atm/s)  in  the  test  chamber.  The  system employs  an  electric 
primer  (FA  874)  as  the  initiator  for  ignition  of  a  propellant  charge  which  is 
used  as  a  source  for  hot  gas  generation.  Using  electric  primers  instead  of  per¬ 
cussion  primers  ensures  higher  reliability  and  greater  reproducibility.  Numerous 
test  firings  have  been  conducted  to  evaluate  the  effect  of  different  types  of 
igniter-propellant  charges  on  the  rate  of  pressurization.  An  ignition  circuit 
has  also  been  designed  and  successfully  used  for  activating  electric  primers. 

4.2  Improvements  in  Data  Acquisition  System 

The  ignition  process  of  solid  propellants  under  highly  transient  pressure 
excursion  is  usually  of  short  duration  (in  the  order  of  a  few  hundred  microseconds). 
To  achieve  better  time  resolution,  a  multichannel  transient  waveform  recorder 
(Physical  Data  Model  #515.234)  was  ordered  in  the  beginning  of  the  contract. 

This  device,  received  in  January  1980,  is  being  used,  together  with  a  magnetic 
tape  recorder  (Hewlett  Packard  Model  3955  B),  to  record  transient  signals.  The 
maximum  sampling  rate  of  the  Physical  Data  System  is  2  x  10^  samples/sec,  and 
the  maximum  amplitude  resolution  is  0.1%.  To  obtain  more  detailed  information 
of  the  ignition  event,  the  *s-frame  optical  head  was  replaced  by  a  %-frame  optical 
head.  Currently,  the  highest  framing  rate  is  about  40,000  pictures  per  second. 

For  extremely  rapid  ignition  events,  streak  photography  was  used.  To  measure 
the  instantaneous  heat  flux  at  the  simulated  propellant  surface,  a  set  of  thin 
film  (0.5  to  1  Urn)  Rhodium-Platinum  heat-flux  gages  (Medtherm  302)  was  ordered 
and  received  only  recently.  These  gages  have  microsecond  response  time  for 
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temperature  measurements.  Due  to  the  delay  in  receiving  the  heat-flux  gages, 
measurements  have  not  yet  been  carried  out  Modifications  are  being  made  in 
the  chamber  in  order  to  install  the  flux  gage. 

4.3  Experimental  Observations  of  Flame  Spreading  and  Combustion  in  Narrow  Cracks 

An  improved  procedure  for  sample  preparations  has  been  developed  to  observe 
and  record  flame-spreading  and  ignition  processes  in  narrow  propellant  cracks. 

A  schematic  diagram  of  the  test  chamber  and  igniter  system  for  this  study  is 
given  in  Figure  8.  A  crack  was  formed  between  a  propellant  slab  and  the  sacri¬ 
ficial  plexiglass  window,  i.e.,  one  side  of  the  crack  was  an  inert,  transparent, 
plexiglass  window,  and  the  other  was  a  propellant  slab  glued  to  a  stainless  steel 
base  plate.  This  type  of  configuration  provided  direct  (front  view)  observation 
of  the  flame-front  propagation  process.  Figure  9  shows  a  typical  propellant 
sample  used  for  detailed  observation  of  the  flame  front.  The  gap  width  of  such 
a  crack  configuration  is  varied  by  the  amount  by  which  the  propellant  surface 
is  recessed  below  the  side-leg  assembly  (see  Figure  9).  The  propellant  slabs 
used  in  these  tests  were  183  mm  long  and  17.7  mm  wide.  Figure  10  gives  a  block 
diagram  of  the  data  acquisition  system. 

Figure  11  presents  a  sequence  of  pictures  of  a  test  performed  on  a  narrow 
crack  (gap  width  *  455  t*n)  in  which  the  crack  was  formed  between  the  propellant 
surface  and  the  plexiglass  window.  It  is  interesting  to  note  that  when  the  width 
of  the  crack  gap  is  very  small  (500  um  or  smaller),  partial  closure  of  the  crack 
occurs  because  of  propellant  deformation.  Partial  closure  first  appears  near 
the  entrance.  As  time  progresses,  the  closure  region  moves  downstream  and  a 
second  partial  closure  region  develops  near  the  entrance.  At  times,  as  the  process 
continues  three  or  four  such  partial  closures  are  observed  simultaneously  (for 
initial  gap  widths  of  the  order  of  450  um).  Usually,  at  the  moment  at  which 


Fig.  10  Block  Diagram  of  Data  Acquisition  System 


Framing  Rate  -  36,700  Pictures  Per  Second 


Fig.  11  A  Sequence  of  Photographs  Showing  Crack-Gap  Closures 
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the  third  closure  develops  at  the  entrance,  the  second  and  first  regions  are 
located,  respectively,  at  approximately  20%  and  50%  of  the  crack  length  from 
the  entrance.  Later,  as  the  combustion  process  becomes  more  pronounced,  the 
partial  closure  regions  disappear  as  a  result  of  both  propellant  regression  and 
higher  pressure  in  the  crack  cavity.  The  entire  process  is  believed  to  be  the 
result  of  deformations  caused  by  high  pressure  acting  on  the  propellant  surface 
exposed  to  the  chamber,  and  by  complex  interaction  between  propellant  deformation 
and  pressure  distribution  in  the  crack  cavity.  This  implies  that  the  coupling 
between  chamber  pressurization,  crack  combustion,  and  propellant  deformation 
is  quite  important,  especially  in  the  case  of  very  narrow  cracks. 

Tests  were  continued  to  study  the  effect  of  pressurization  rate,  propellant 
type,  crack-gap  width,  and  crack  length  on  flame-spreading  and  combustion  processes 
in  solid  propellant  cracks.  Details  of  the  procedure  and  the  results  are  given 
in  Ref.  18.  The  abstract  for  the  paper  is  given  in  Appendix  IV. 
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V.  NOMENCLATURE 


A 


c>  cp 


D 

E 


G 


r 


t 

v 

z 


z 


frequency  factor 
specific  heat 

binary  mass  diffusion  coefficient 
activation  energy 
black  body  radiation  intensity 
external  radiation  heat  flux 
local  radiation  flux  at  z 

thermal  conductivity 

rate  constants  for  reactions  1  and  2 

pressure 

heat  flux  (energy  per  unit  time  per  unit  area) 

heat  generation  rate  (energy  per  unit  time  per  unit  volume) 

rate  of  heat  generation  due  to  photochemical  process  per 
unit  wavelength/local  radiation  flux  at  z 

heat  of  pyrolysis  per  unit  mass 

radial  distance  from  the  center  of  the  statistically 
averaged  element 

burning  rate 

gas  constant 

universal  gas  constant 

radii  of  outer  surfaces  of  oxidizer  particle  and  fuel 
binder,  respectively 

time 

gas  phase  velocity  in  z  direction 

mass  fraction  of  species  i  in  the  gas  phase 

distance  from  the  initial  (t*0)  position  of  the  interface; 
positive  in  the  gas  phase 


a 

6 


nX 

X 

P 

T 


radiation  absorptivity  of  the  interface 
in-depth  radiation  absorption  coefficient 
emissivity  of  the  interface 

[~ 2  2 

radial  distance  of  the  spherical  coordinate,  r  +  z 
efficiency  for  photochemical  process  at  wavelength  \ 
wavelength  of  external  radiation 
density 

transmissivity  of  the  solid  phase 
mass  production  rate  per  unit  area 
mass  production  rate  per  unit  volume 


Subscripts 
1 ,2,a,b,c 
F 

F-g 

g 

i 

Ox 

Ox-F 

Ox-g 

pc 


py 

s 

z 


X 


various  reactions 
fuel 

fuel-gas  interface 
gas  phase 

initial  value  or  species  i 
oxidizer 

oxidizer-fuel  interface 
oxidizer-gas  interface 
photochemical 
pyrolysis 
solid  phase 

in  z  direction  (perpendicular  to  the  initial  position 
of  the  interface) 

far  away  from  the  interface  in  z  direction 
wavelength 
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APPENDIX  I 


REVIEW  OF  SOLID-PROPELLANT  IGNITION  STUDIES 


A.  K.  Kulkaml,  M.  Kumar ,  and  K.  K.  Kuo 
Department  of  Mechanical  Engineering 
The  Pennsylvania  State  University 
University  Park,  PA  16802 


ABSTRACT  (Ref.  2) 

An  extensive  review  of  the  literature  on  solid-propellant  ignition 
has  been  made  to  establish  the  state  of  the  art.  Various  ignition  theories, 
experimental  measurements,  and  ignition  criteria  are  critically  examined. 

The  review  is  summarized  In  an  easy-to-read  tabular  form  to  facilitate 
comparison  between  various  studies.  The  effects  of  important  parameters 
on  ignition  processes  are  also  discussed.  Major  technological  gaps  are 
identified,  and  areas  for  future  studies  are  recommended. 
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APPENDIX  II 

COMBUSTION  -  STRUCTURAL  INTERACTION  IN  A  VISCOELASTIC  MATERIAL 

T.  Y.  Chang  and  J.  P.  Chang 
Department  of  Civil  Engineering 
The  University  of  Akron,  Akron,  Ohio 

M.  Kumar  and  K.  K.  Kuo 
Department  of  Mechanical  Engineering 
The  Pennsylvania  State  University 
University  Park,  Pennsylvania 


ABSTRACT  (Ref.  16) 

The  effect  of  interaction  between  combustion  processes  and  structural 
deformation  of  solid  propellant  was  considered.  The  combustion  analysis 
was  performed  on  the  basis  of  deformed  crack  geometry,  which  was  deter¬ 
mined  from  the  structural  analysis.  On  the  other  hand,  input  data  for  the 
structural  analysis,  such  as  pressure  distribution  along  the  crack  boundary 
and  ablation  velocity  of  the  crack,  were  determined  from  the  combustion 
analysis.  The  interaction  analysis  was  conducted  by  combining  two  computer 
codes,  a  combustion  analysis  code  and  a  general  purpose  finite  element 
structural  analysis  code. 
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APPENDIX  III 


EFFECT  OF  PROPELLANT  DEFORMATION  ON  IGNITION 
AND  COMBUSTION  PROCESSES  IN  SOLID  PROPELLANT  CRACKS 


M.  Kumar  and  K.  K.  Kuo 
Dept,  of  Mechanical  Engineering 
The  Pennsylvania  State  University 
University  Park,  PA  16802 


To  be  Presented  at  The  1980  JANNAF 
Propulsion  Systems  Hazards  Meeting 
to  be  held  at  the  Naval  Postgraduate 
School,  Monterey,  California,  October 
27-28,  1980 
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ABSTRACT  (Ref.  17) 


In  order  to  obtain  a  detailed  observation  of  the  ignition-front  propagation 
in  narrow  solid  propellant  cracks,  tests  were  performed  on  cracks  formed  between 
a  propellant  slab  and  an  inert,  transparent  plexiglass  window.  Tests  were  con¬ 
ducted  under  rapidly  increasing  chamber  pressure.  Ignition-front  propagation 
rates  were  measured  using  high-speed  (up  to  44,000  pictures  per  second)  camera. 

It  was  observed  that  when  the  width  of  the  crack  gap  is  very  small  (<500  ym) , 
partial  closure  of  the  crack  occurs  due  to  propellant  deformation.  This  partial 
closure  first  appears  near  the  entrance.  As  time  progresses,  the  closure  region 
moves  downstream  and  a  second  partial  closure  region  develops  near  the  entrance. 

This  process  continues;  at  times,  three  or  four  such  partial  closures  are  observed 
simultaneously  (for  initial  gap  widths  of  the  order  of  450  ym) .  Usually,  at  the 
moment  when  the  third  closure  develops  at  the  entrance,  the  second  and  first  are 
located,  respectively,  at  approximately  20%  and  50%  of  the  crack  length  from  the 
entrance.  Later  as  the  combustion  process  becomes  more  pronounced,  the  partial 
closure  regions  disappear  as  a  result  of  both  propellant  regression  and  higher 
pressure  in  the  cavity. 

The  entire  process  is  believed  to  be  the  result  of  deformation  caused  by 
high  pressure  acting  on  the  propellant  surface  exposed  to  the  chamber,  and  by 
the  complex  interaction  between  propellant  deformation  and  pressure  distribution 
in  the  crack  cavity.  Propellant  surface  exposed  to  the  chamber  are  compressed 
by  the  high  chamber  pressure,  which  results  in  the  propellant  being  pushed  into 
the  crack.  Since  during  the  initial  pressure  transient,  the  chamber  pressure 
increases  faster  than  the  pressure  in  the  crack,  the  propellant  is  pushed  toward 
the  lower  pressure  region  inside  the  crack.  The  mechanical  deformation  of  the 
propellant  causes  narrowing  of  the  crack  gap  and  consequently  results  in  local 
crack  closure. 

The  interaction  between  structural  deformation  a'nd  combustion  processes  was 
also  investigated  theoretically.  In  the  theoretical  model,  the  effect  of  inter¬ 
related  structural  deformation  and  combustion  phenomena  was  taken  into  account 
by  considering:  a)  transient  one-dimensional  mass,  momentum  and  energy  conserva¬ 
tion  equations  in  the  gas  phase,  b)  a  transient  one-dimensional  heat  conduction 
equation  in  the  solid  phase,  and  c)  quasi-static  deformation  of  the  two-dimensional, 
linear- viscoelastic  propellant  crack  due  to  pressure  loading.  This  set  of  coupled, 
nonlinear,  partial  differential  equations  was  solved  numerically.  The  gas-dynamic 
equations  were  solved  using  a  finite  difference  analysis  and  the  structure  mechanics 
equation  was  solved  using  a  finite  element  analysis.  Calculated  results  of  pressure 
distributions  at  various  times  indicate  that  the  effect  of  propellant  deformation 
is  quite  important.  The  gap  closure  results  in  substantial  local  pressure  peaks. 

Results  obtained  Indicate  that  closure  of  the  crack  gap  initially  occurs 
at  the  crack  entrance  because  of  small  gap  widths  and  high  chamber  pressures. 

The  gap  closure  is  not  observed  for  cracks  with  large  gap  widths  or  at  low  chamber 
pressure,  because  propellant  deformation  under  these  conditions  is  small  and  does 
not  affect  the  ignition  process  substantially.  Gap  closures  and  resulting  pressure 
peaks  can  strongly  influence  the  convective  burning  process  in  the  crack  and  may 
contribute  to  the  deflagration-to-detonation  transition  process. 
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APPENDIX  IV 


GAS  PENETRATION,  FLAME  PROPAGATION,  AND  COMBUSTION  PROCESSES 
IN  SOLID  PROPELLANT  CRACKS 

Mridul  Kumar,  Stephen  M.  Kovacic,  and  Kenneth  K.  Kuo 
The  Pennsylvania  State  University 
University  Park,  Pennsylvania 


ABSTRACT  (Ref.  18) 

The  effects  of  pressurization  rate,  crack-gap  width,  crack  length, 
and  propellant  type  on  the  ignition  and  flame-spreading  processes  in 
isolated  AP-based  solid  propellant  cracks  have  been  studied  experimen¬ 
tally.  Ignition  front  propagation  rates  were  measured  using  a  high¬ 
speed  (up  to  44,000  pictures/second)  camera.  Cracks  up  to  200  mm  in 
length  with  gap  widths  as  low  as  450  pm  were  studied.  It  was  observed 
that  the  hot  gases  precede  the  ignition  front.  The  ignition-front  pro¬ 
pagation  speed  increases  near  the  crack  entrance,  reaches  a  maximum, 
and  then  decreases  near  the  crack  tip.  The  results  of  parametric  study 
indicate  that  the  time  required  for  the  ignition  front  to  reach  the 
crack  tip  decreases,  and  that  the  maximum  velocity  of  the  ignition  front 
increases  as  the  pressurization  rate  or  burning  rate  of  the  propellant 
is  increased.  The  maximum  pressure  in  the  crack  increases  with  an  in¬ 
crease  in  burning  rate  or  crack  length,  but  decreases  with  an  increase 
in  gap  width.  It  was  observed  that  when  the  crack  gap  width  is  very 
small  (~450um) ,  several  partial  closures  of  the  crack  occur  due  to  pro¬ 
pellant  deformation  caused  by  the  strong  interaction  between  combustion 
and  structure  mechanics. 
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